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Abstract 

We have derived several relations, which allow the evaluation of the system free energy changes in 
the leading order in h"^ along classically generated trajectories. The results are formulated in terms 
of purely classical Hamiltonians and trajectories, so that semiclassical partition functions can be 
computed, e.g., via classical molecular dynamics simulations. The Hamiltonians, however, contain 
additional potential-energy terms, which are proportional to and are temperature-dependent. 
We discussed the influence of quantum interference on the nonequilibrium work and problems with 
unambiguous definition of the semiclassical work operator. 



1 



I. INTRODUCTION 



The nonequilibrium work theorem, or Jarzynski equation, relates nonequihbrium work 
performed on a dissipative system during certain time interval to the difference of the corre- 

the work 
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theorem has been generalized to the quantum case, too 
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Monnai and Tasaki have found quantum (~ fi^) corrections to 



the transient and steady-state fluctuation theorems for a damped harmonic oscillator [21 1. 
Chernyak and Mukamel gave quantum (~ h^) corrections to the work theorem for a spe- 
cific form of the driven system Hamiltonian 22]. The quantum generalizations of the work 
theorem have triggered a certain controversy in the literature on whether it is possible to 
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uel y d efine the quantum work operator and how to interpret the quantum work theorem 
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In the present paper, we use the ideas inspired by the classical work theorems to de- 
rive several relations which allow us the calculation of equilibrium semiclassical free energy 
changes along non-equilibrium classical trajectories. Our aim is twofold. First, we develop 
a practical tool for the calculation of semiclassical free energies differences through classical 
equations of motion and/or molecular dynamics simulations. Such semiclassical corrections 
are responsible, e.g., for the short-distance anomaly of the pair distribution function of liquid 



neon 



23 



241 ]. Second, we wish to get a better understanding of how quantum interference 



influences the nonequilibrium work and if it is possible to uniquely define the semiclassical 
work operator. 

Note that all the semiclassical quantities which are considered in the present paper are 
evaluated in the leading order in h"^, i.e., with the accuracy 0{h'^). We use the simbols h 
and % intermittently throughout the text, in order to avoid writing numerous factors of 27r. 

II. WIGNER DISTRIBUTIONS AND USEFUL IDENTITIES 



Let us consider a collection of quantum point particles with the Hamiltonian 



N ^2 
1=1 



(1) 
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Qii Pi = —ihd/dqi, and rrii are the positions, momenta, and masses of the particles, U{q) is 
the potential energy. Hereafter, the boldface notation p, q is used to collectively denote the 
set of all Pi and gj. We assume that the particles are prepared in an equilibrium canonical 
ensemble at the temperature T. The corresponding distribution (density matrix) reads 

/5(p, q) = exp{-/3i^(p, q)}, Z = Tr(exp{-/?i7(p, q)}), (2) 

Z being the partition function. 

We wish to calculate the leading (~ h^) corrections to the canonical distribution ([2]) and 
partition function Z. To this end, it is convenient to switch to the Wigner representation 



25|, l26|, so that ([2j) becomes the corresponding Wigner distribution (the subscript W) 

p[J)(p, q) = Z-' exp{-P{H{p, q) + h'A^'\p, p, q))} + 0{h'), (3) 

Z = J dprfqexp{-/3(if(p,q) +/i2AW(/3,p,q))}. (4) 

Here p and q can be treated as the phase variables in the Wigner space, taking the trace 
reduces to the integration over the Wigner phase space variables, -^(p, q) is the classical 
Hamiltonian corresponding to its quantum counterpart ([T]), and A'^^-' is the temperature- 



dependent quantum correction 
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For our further purposes we preaverage over momenta. That is, we replace PiPj by 
^ijirii/ (3 and obtain 



i=l 

We also introduce two new functions 

N 



4f 1 12m, dq} 24mi V % / I ' 



Z^\24m, dql J' 



i=l 

The physical meaning of these new functions will be explained below. Let us now consider 
the quantities 

p[^)(p, q) = Z-' exp{-P{H{p, q) + h'A^^\p, p, q))}, a = 0, 1, 2, 3. (9) 



Here is the semiclassical Wigner distribution ([3]) but p[^^ for a = 1, 2, 3 are not true semi- 
classical Wigner distributions [27,]. However, all p^j share the following evident property: 
The partition functions 

Z(») = y"rfprfqexp{-/3(i7(p,q) + /i2A(»)(AP,q))} (10) 

coincide with the true semiclassical partition function Z^^"* within the accuracy Oiji^). The 
equivalence of Z^^-* and Z^'^^ , Z^^^ can be demonstrated via integration by parts, 



P 



J dqexp{-f3U{^)} (^^) = / dqexp{-(3U{^)}- 



dql 



III. SEMICLASSICAL FREE ENERGY CHANGES ALONG CLASSICAL TRA- 
JECTORIES 

We are in a position now to derive the semiclassical version of the nonequilibrium work 
theorem. Let 

N „2 

^(p,q,t) = 5^^ + f/(q,t). (11) 

be a quantum Hamiltonian, which is allowed to be explicitly time-dependent. Let us now 
introduce the quantities 

H^'^^ p, q, t) = H{p, q, t) + /i^ A^") p, q, t) , a = 0, 1, 2, 3. (12) 

Here H{p, q, t) is the classical Hamiltonian which corresponds to the quantum Hamiltonian 
( HT]) . and functions A^^\P, p, q, t) are constructed out of H{p, q, t) as is prescribed via Eqs. 

m-m- 

Let us now consider H^"''^ for a = 1, 2, 3 as classical Hamiltonians which are parametrically 
temperature-dependent j28|. They can equivalently be rewritten as 

N 2 

H^^^ p, q, t) = 5^ + {f3, q, t) , f/^'^) {(3, q, t) = U{q, t) + h'^A^'^^ q, t) . (13) 
^-^ Irrii 

i=l 

So, the quantum correction gives rise to a simple redefinition of the potential energy, which 
acquires an additional temperature-dependent ~ fi^ contribution. Now we can use the 
Hamiltonian (fT3l) to write down the corresponding classical equations of motion 

dp,{t) _ gW(/?,p,q,t) dq,{t) _ gW(/?,p,q,t) 
dt dqi ' dt dpi 



which produce the classical trajectories p(t), q(i)- Furthermore, we can treat 

pi^)(p,q,t) = Zr^exp{-/?i7('^)(/?,p,q,t)} (15) 

as the corresponding classical canonical distribution with the partition function 

Zt = J dpdciexp{-(3H^''\(3,p,q,t)}. (16) 

After this is done, we can proceed analogously to the derivation of the classical Jarzynski 
equation 0, Q . Namely, we can write the following chain of identities: 

dpodqoexp{-PH^'^\P, po, qo, 0)} exp{-P{H<^^\p, p^, q^, t) - H^-\f5, po, qo, 0))} = 

j dpodqoexp{-l3H^''\(3,pt,qut)} = j dpirfq* exp{-/5i/('^)(/5, p^, q^, t)} = Zt, (17) 

a = 1,2, 3. In deriving Eq. ( |T71) we make use of the fact that motion of a Hamiltonian 
system can be regarded as the canonical transformation, for which the Liouville theorem 
holds: dpodqQ = dptdqt- Introducing the abbreviation 

J dpodqoZ^' exp{-(3H^^\(3,po,qo,0)}.: = (••■)o (18) 
and dividing Eq. f|T7j) by Zq, we obtain 

(exp{-/?(/7('^)(/3,pt,qt,t) - ff(")(APo,qo,0))}>o = Z^/Zo, a = 1,2,3. (19) 

We can further split the total Hamiltonian into the system Hamiltonian, the bath Hamil- 
tonian and their coupling, 

H^''\P, p, q, t) = H^iP, PS, qs, t) + Hf{(3, p^, qs) + Ps> ^s, Pb, q^)- (20) 

Here the system Hamiltonian only is allowed to be explicitly time-dependent. Plugging H'^'^^ 
( |20|) into the identity (fT9!) and making use of the equations of motion (|T^ . we can write 

//(")(/?, Pi, qi,t)-//('^) (A Po,qo,0) = 

iyt'-^H^^\P,p{t'),qit'),t') = j\t'^Hf{P,ps{t'),qs{t'),t') = W^^\ (21) 
W^""^ being the work performed on our classical system. Thus Eqs. ( |T9l) and ( 12T1) yield 

(exp(-/3W^('^))>Q = Z,/Z,, (22) 
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which is the semiclassical analogue of the classical Jarzynski formula 

Note that the ratio of the total {S + B) partition functions Zt/Zo equals to the ratio of 
the system partition functions k Meed, the reduced system (S) density matrixiiis 



defined via averaging the total density matrix f[T51) over the bath degrees of freedom 29|, |30| : 



pP{Ps,^s,t) = Zsl j rfpBC^qsexp{-/3(i7('^)(/?,p5,qs,PB,qB) -i^ij)}, (23) 

where the bath free energy Fb = — {In Zb)/ (3 is determined through the bath partition 
function 

Zb = j rfpBc/qijexp{-/?if^"^(/?,Pij,qB)}. (24) 

The system partition function in Eq. (!23l) is explicitly defined as Zs^t = Zi/Zb and we can 
write 

Zt/Zo = Zs,t/Zs,o- (25) 



If the quantum ~ h corrections are neglected, then Eos. (1191) and ( ]22l) reduce, of course, 
to the classical work theorem for Hamiltonian systems 3|. 

To illustrate computational aspects of practical use of relations fll9|22p . we performed 
"molecular dynamics" simulations for a point particle in one-dimensional nonlinear potential. 
The quantum Hamiltonian reads 

H{p,q,t) = ^ + af + bq'j^^. (26) 

In this example, the mass m = 1, the parameters a = b = 1, the position q, the momentum 
p = —ihd/dq, the time t, and the Planck constant are taken as dimensionless. We con- 
struct the corresponding semiclassical Hamiltonian H^^\P,p, q, t) as is explained in Sections 
II and III. Then we calculate the free energy difference at times t and 0, AF = Ft — Fq = 
— \n.{Zt/Zo)/f3, by "molecular dynamics" simulations and "exactly". In "molecular dynam- 
ics" simulations, we sample initial values of q and p according to the Wigner distribution 
(fT5|) {Nsam being the number of samplings), solve equations of motions (IT^ numerically, and 
calculate the quantity exp{— /?(if'^^)(/?,p, g, t) — H^^\l3,p, q, 0))}. The procedure is repeated 
Nsam times, the averag ed value {exp{-p{H^^\p,p,q,t) - H'-^\p,p,q,0))})^ is obtained, 
and AF is finally calculated through the semiclassical work theorem f[T^ . "Exact" AF is 
obtained trough the direct numerical evaluation of the semiclassial partition functions Zt 
and Zq according to Eq. (fT6|) . The results of the calculations are depicted in figure 1. We 
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see that making Ngam, = 10^ samplings gives already reasonable, but quite noisy estimation 
for AF (dashed lines), Ngam = 10^ gets a better result (dotted lines). Nsam = 10^ yields 
AF which is virtually indistinguishable with the exact results (full lines). The convergence 
is thus rather slow (see Ref. {si]] for the comparison of different simulation schemes for 
obtaining free energy differences). On the other hand, the convergence speed is the same 
for classical (bottom curves, h = 0) and semiclassical (upper curves, h = 1) case 32|. 

The Hamiltonian H^^\P,p,q), a = 0, has already been used in molecular dynamics 



simulations 



33|. Such a choice gives a correct semiclassical distribution, but its practical 



implementation may encounter certain difficulties 28|, |33|]. The choice a = 1 gives a correct 
semiclassical distribution over positions q. The corresponding Hamiltonian H^^\P, p, q) has 
also found its application in molecular dynamics simulations, being evaluated up to the terms 
of the order of |23| and |2J, ^ . To our knowledge, the Hamiltonians H^'^\P,p,q) 



and H^^\P,p,q) have never been used in practice. However, they give more tractable 
expressions for the semiclassical potential A*^"''(q, t). It should be noted that (almost) all 
physically relevant potentials U{q) when |q| —>■ oo. In that case, there is no problem 
with the exponentiating the semiclassical correction potentials, e.g.. 



exp{-l3H{p, q, t)}{l + h'(3A^'^^ (/3, p, q, t)) = e^-f^H'^KP: P, ^^t)) + 0{h'). (27) 

Such exponentiation may work even beyond its strict domain of validity 0(^^), since it 
corresponds to the partial summation of the higher order contributions in the semiclassical 
expansion. If we use strongly attractive potentials, then exponentiation in Eq. (1271) for 
a = 1 may not be feasible, since the contribution due to the second term in Eq. ([6]) may 
become predominant, so that the exponential explodes. In such a case, the choices a = 2 
and a = 3 are necessary. Finally, a word of caution concerning applications of the present 
(and similar) semiclassical methods to realistic molecular systems. The ubiquitous presence 
of sharply varying repulsions in condensed phases means that the terms like ([7j) and ([8]) can 
be quite large. This can lead to numerical difficulties jssl. More fundamentally, this could 
indicate that a straightforward ^^-perturbative treatment might not always be appropriate. 



7 



IV. DIFFICULTIES WITH THE DEFINITION OF SEMICLASSICAL WORK 



Several important points are to be discussed here in connection with the interpretation 
of our results as semiclassical work theorems. As has been stressed in Sec. II, the distribu- 
tions p|^''(p, q) ([2]) for a = 1,2,3 are not true semiclassical Wigner distributions. However, 
they give correct values of the semiclassical partition functions within the accuracy 0{h'^). 
Furthermore, the Hamiltonian dynamics governed by Eqs. f|T4|) is not, of course, the true 
semiclassical dynamics. However, if we run classical molecular dynamics simulations with 
the Hamiltonians p, q, t) (a = 1,2,3) we get correct semiclassical values of the par- 

tition functions, again with the accuracy 0{h^). 

The explicit form of the semiclassical nonequilibrium "work theorem" (Eqs. (HM and 
f l22|) ) is not unique. We have simultaneously derived four (a = 0,1,2,3) different "work 
theorems" with different definitions of the work operator W^""^ (1211) . For each particular 
choice of the classical Hamiltonian, W^""^ can definitely be treated as the classical work, 
and all W^""^ coincide in the limit h ^ 0. However, we cannot regard any of W^""^ as a 
true semiclassical work operator. First, W^""^ for a = 0, 1, 2, 3 yield different values of work, 
and it is not clear which choice (if any) is preferable. Second, the mean value of the work 
operator W^"''' does not coincide with the mean semiclassical energy difference AE. Indeed, 
adopting the notation 



j dptdqtZi'exp{-PH^^\p,pt,c^t,t)}... = (...),, 



(28) 



we can explicitly define 



(29) 



and 



AE = {HiP, pt, qt, t)), - {HiP, po, qo, 0))^ . 



(30) 



Therefore, 



W'^''^ = AE + h' {(A('^)(Aqt,t)>, - (A('^)(Aqo,0)> J . 



(31) 



The difference between W^""^ and AE is of the order of h"^ and is generally nonzero 35|]. This 



is the semiclassical way of stating t 



define the quantum work operator [10|: 



lat it is hardly possible to properly and unambiguously 

iil,[i2l,Q,Q,0,Q. 
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It is important that W^"''' fl^T]) is explicitly temperature dependent. That is a direct con- 
sequence of quantum interference. In classical case, we can consider a single trajectory in 
the phase space and calculate the work along it. Such a work is determined exclusively by 
the difference of energies along the trajectory. Incoherent summation over many such tra- 
jectories, weighted via classical canonical distribution, yields the classical work theorem. In 
quantum or even semiclassical case, there exists interference between the trajectories, which 
manifests itself through the temperature dependent Hamiltonians H^"-\P,p,q,t). Thus, 
even in the semiclassical limit, it is not possible to define W^""^ as a function of exclusively 
the phase space variables and time. The ensemble quantity, viz., the temperature, inevitably 
enters the definition of work. 

The semiclassical Jarzynski equations f|T9l) and fl22|) cannot be directly considered as 



semiclassical versions of the quantum Jarzynski equation 



14 



16 



17 



, Q, Q, Q . Indeed, the latter is written as |ll[ | 
exp{l3H{0))exp{-l3H{t)) 



iQ, y, y, 0, y, s 



10 



11 



12 



13|, 



T<(exp(-/5#)) = V^o- 



(32) 



Here H{t) is the total Hamiltonian in the Heisenberg representation, (...) 



Tr(exp{— /3if(0)}...), and Ty is the chronologica 



dependent work operator is explicitly defined as 



d 



13, 



ordering operator. The trajectory- 



18| 



W = H{t) - HiO) = I dt'—Hit') = I dt'—Hsit') 



d 



dt' 



(33) 



(the last expression in Eq. ([3SD holds true provided we partition the total Hamiltonian as a 
sum of the system Hamiltonian, bath Hamiltonian and their coupling, and assume that the 
system Hamiltonian is explicitly time dependent, H{t) = Hs(t) + Hb + Hsb)- In principle, 
we can write down a direct semiclassical analogue of Eq. fl32l) . The corresponding formulas 
are not presented here, since they are quite cumbersome and difficult to interpret. What is 
even more important, we have to introduce the semiclassical time-evolution operator, which 
should be evaluated in the leading order in . That would prevent us from a simple classical 
interpretation of the formulas. 



V. CONCLUSIONS 



We have derived semiclassical analogues (fT9l) and (p2l) of the classical work theorems, 
which allow us to evaluate the ratio of the (system) partition functions in the leading order 

9 



in h'^, with the accuracy 0{h^). The corrections due to the quantum statistics (Bose-Einstein 
or Fermi-Dirac) give rise to the contributions ~ to the equihbrium density matrices and 
partition functions jsG] and are not considered here. 

The semiclassical analogues of the work theorems are formulated in terms of purely classi- 
cal Hamiltonians if^"-* (/?, p, q) (fT2|) and classical trajectories, so that semiclassical partition 
functions can be evaluated, e.g., via classical molecular dynamics simulations. The Hamilto- 
nians, however, contain additional potential-energy terms, which are proportional to fi^ and 
are temperature-dependent. The Hamiltonians H^'^\P,p,q) and H^^\P,p,q) have already 
been used in classical, the so-called Wigner-Kirkwood, molecular dynamics simulations (see 



Refs. 



33|, l37| and 



3, 



24 



341, correspondingly). Conceptually similar are also the Feynman- 



Hibbs molecular dynamics j38| and the quantized Hamilton dynamics 39|. We mention also 
that the use of the path integral technique in dynamical computer simulations necessitates 
introducing additional fictitious classical degrees of freedom into the system under study. 



For example, the riiig 



polymer beads 



37 



jolymer molecular dynamics is based on the classical evolution of the 



40 



4l| . The path integrals have recently beei i ap plied to the calcula- 



431] . In our semiclassical 



tion of partition functions through the quantum work theorem |42l . 
method, no additional degrees of freedom are introduced. Furthermore, the higher order 
in contributions can straightforwardly be incorporated into our main formulas ( fT9|) and 



221) . if necessary. For examp 



has been used in Refs. 



24 



34|. 



e, p, q) evaluated up to the terms of the order of 



The explicit form of the semiclassical work theorem (Eqs. f|T9|) and fl22|) ) is not unique. In 
fact, we have derived four (a = 0,1,2,3) different work theorems with different definitions 
of the work operator W^'^^ fl2T|) . For each particular choice of the classical Hamiltonian 
H^"'\P,p,q), W^""^ can definitely be treated as the classical work, and all W^'^^ coincide in 
the limit h 0. However, we cannot claim any of W^""^ to be a true semiclassical work 
operator. First, W^""^ yield different values of the work performed, and it is not clear which 
choice (if any) is preferable. Second, the mean value of the work operator W^^^ (1291) does not 
generally coincide with the mean semiclassical energy difference AE fl30l) . This indicates, at 
the semiclassical level, that it i s hardly possible to properly and unambiguously define the 



3, 



11 



12 



13 



14 



13, 



18(1 . However, the ambiguity in the physical 



quantum work operator 
meaning of the quantity W^""^ does not prevent us from using the semiclassical work theorem 
in practical calculations and classical molecular dynamics simulations, in order to get the 
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system free energy differences or semiclassical averages. 

Note finally that the present approach can be incorporated into a general scheme devel- 
oped in to generate semiclassical versions of different fluctuation theorems for Hamilto- 
nian systems in the equilibrium and in a steady state. For example, treating the Hamilto- 
nians fl20p as true classical Hamiltonians, we can straightforwardly derive the semiclassical 
counterpart of the Crooks transient fluctuation theorem 

{6{w - lV('^nPo,qo,t))>o = + W^^\pt,c[t,t)))^exp{(3w)Zt/Zo. (34) 

Here W^"'\po,<lo,'t) and W^'^\pt,(lt,'t) is the classical work (pTl) expressed as the function 
of the initial and final coordinates, correspondingly. 
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0.25 n 




t, reduced units 

FIG. 1: The difference between the free energies at times t and 0, AF = Ft — Fq = — ln(Zt/Zo) / /3, 
as a function of t. The lower curves show the classical result {h = 0) and the upper curves show 
the semiclassical result {h = 1). Full lines correspond to "exact" semiclassical calculations, dashed 
and dotted lines correspond to "molecular dynamics" simulations (see text for details). 
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